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ABSTRACT 

The motion of an unbalanced rotor during acceleration through its 
critical speed is studied by numerical solutions obtained with a digi- 
tal computer. The rotor is laterally restrained in two orthogonal 
directions by linear springs and is accelerated by a constant applied 
torque. It is found that, for a fixed combination of rotor unbalance 
and lateral stiffnesses, the applied torque must exceed a limiting 
value in order to accelerate successfully through the critical speed 
region. With smaller applied torques, after initial acceleration from 
rest the speed oscillates continually about the critical speed and the 
lateral excursions grow steadily. The conditions necessary for a 
successful acceleration are established and the maximum lateral ex- 
cursions during successful accelerations are determined. The effects 
of small amounts of viscous damping in the lateral directions are also 
obtained. Finally, an example problem is solved to illustrate the 
extension of the solution to those problems where applied torque is 


a function of speed. 


il 








ACKNOWLEDGMENT 
The writer wishes to express his appreciation to Dr, Robert E. 
Newton for providing the stimulus which initiated this investigation, 
and whose guidance and encouragement aided immeasurably in its success- 


ful completion. 


iiL 





Section 


Abstract 


Acknowledgment 


Table of Contents 


Notation 


Bibliography 
Illustrations 


Appendix I 


Appendix II 
Appendix III 
Appendix IV 


Appendix V 


TABLE OF CONTENTS 


Title 


Introduction 

Equations of Motion 

Numerical Solution 

Computer Program -Design and Use 
Discussion of Results 


Sample Problem 


Conversion of Equations for Numerical 
Integration 


Selection of Time Increment 
Block Diagram 
Computer Frogram 


Tabulated Data 


iv 


Page 
ii 
iii 


iv 


10 
15 
20 
22 
23 


43 


47 
48 
51 


54 





P(1) 
P(2) 
P(3) 
P(^) 
267 


R/e 


x,y 


09% 


1 


" 


H 


0 0086 
Damping in x direction 
Critical damping in x direction 
Damping in y direction 
Critical damping in y direction 
Eccentricity of mass center from center of rotation 


Moment of inertia about a longitudinal axis through the center 
of rotation 


Moment of inertia about a longitudinal axis through the mass 
center 


Combined shaft and bearing stiffness in the x direction 
Combined shaft and bearing stiffness in the y direction 
Pounds force 

Pounds mass 

Mass of the system 

Applied torque 

Eccentricity parameter mehr 

Torque parameter Mm/Ik, 

Stiffness ratio k,/k, 

Damping parameter (2c,/c رت‎ = (2c,/c , (ى‎ 

Modifying parameter for applied torque 

Dimensionless resultant deflection 

Dimensionless deflection in the x direction 
Dimensionless deflection in the y direction 

Angle measured counter clockwise from x axis to radial axis 


originating at the center of rotation and passing through the 
mass center 


Dimensionless time 
Coordinates of center of rotation 


Coordinates of mass center 





1]. Introduction. 

The maximum amplitude of vibration of an unbalanced rotor upon 
acceleration through its critical speed is of real importance to design 
engineers. 

Àn analytic solution was obtained by Lewis 11 for the linear 
single degree of freedom system having constant acceleration and con- 
stant force amplitude, with and without damping. Baker [2] obtained a 
mathematical machine solution of a linear two degree of freedom system 
with constant acceleration, also with and without damping. Meuser and 
Weibel [3] obtained a solution on the mechanical analyzer for a single 
degree of freedom system having constant acceleration and linear plus 
cubic elasticity, with and without damping. An analog computer solution 
to the damped and undamped linear two degree of freedom system was ob- 
tained by McCann and Bennett [4] . Here again a constant acceleration 
was assumed, Dornig [5] solved analytically the undamped single degree 
of freedom system with constant acceleration. 

Both Biezeno and Grammel n and Baker [2] mention the possibility 
of an unsuccessful acceleration; that is, an acceleration to the vicinity 
of but not through the critical speed. In this case, the energy supplied 
is absorbed by the vibrations and damping, if present, rather than be- 
ing absorbed by the rotor itself in the form of kinetic energy of ro- 
tation. 

In the material that follows, a system having two lateral degrees 
of freedom is investigated. The rotor is accelerated by a constant 
applied torque rather than having a constant angular acceleration. For 
this system, the maximum amplitudes of vibration will be obtained for an 
appropriate range of dimensionless parameters. Included will be the 
damped and undamped cases, and those of equal and unequal spring con- 


stants in the x and y directions. In addition, the areas of successful 


* Square brackets refer to bibliography 








and unsuccessful acceleration will be defined, Finally, an example 
will illustrate the extension of the solution to those problems where 


torque is a function of speed. 








2. Equations of Motion. 

Figure 1 “shows the idealized system, The rotating mass m has 
moment of inertia I, about the longitudinal axis through the mass 
center which is displaced a distance e from the center of the support- 
ing shaft S. Lateral displacement is resisted by the stiffnesses of 
the shaft and its bearings, k, and k, and the associated viscous damp- 
ing cy and ei A driving torque M accelerates the rotor. 

Figure 2 shows the free body diagram of the system and includes 
all real and D'Alembert forces. The origin of the coordinate system 
is at point B, which is the location of the shaft center S when the 
rotor is undeflected. 

Summing forces and moments: 

A’ d 


SE سے‎ = T Rix = mel) 


Er at = 





POA Bus eb 72) همه‎ O (2) 


| ae d 4B v? 
ZFy: iis + C2 +- R2 4y E AL? 


At) 


dt 2 2 


© سر رت ہہ 


EUM FE = — IY و۳‎ 
D S Le M € tdt? 


Taking the moment of inertia I about the shaft center S rather 


than the mass center reduces Eq. 3 to: 


1*6 de E 
o mel — 4 
Lom =M- nel TTE) Ca B + mel i) di 6 (4) 


3 
where [= |, + me”, 
It is convenient at this point to reduce Eqs. 1 and 2 to dimension- 


less form by dividing by ek,, and Eq. 4 to dimensionless form by divid- 


l 
۱ 2 

ing by e ki. 
Then Eq. 1 becomes: 


*Figures begin on page 23. 








A A mide d'en ۰ 
na (4) ex s m LI + log (5) 


Eq. 2 becomes: 


(P ouo (5%) م‎ | (6) 
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and Eq. 4 becomes: 


I E e EE At, IM [ls [9] ا‎ (7) 











ety dt? e E ek, ALT خرل‎ 
Let = dimensionless deflection in the x direction = = 
N = dimensionless deflection in the y direction = + 
N 0۶ 
Ë = dimensionless time = nk 
dS 77٤6 b, 
ii om = f — where the dot denotes differentia- 
"ns AS? dt m 7 
tion with respect to dimensionless time 5 e 
Then Eq. 5 becomes: 
> E E 2512. 27 ` 
s y +J = (8) سم ۵ + 6 میع‎ ۵ (8) 
N nk 
Eq. 6 becomes: 
oe Ca : .ı2 ' ee 
vane 4 + ty =@) وعم‎ - 9 Cod O (9) 
N m kı £, 
me* 
and Eq. 7 becomes, after multiplication by : 
et Mm e? SR ۱ oe 
ee f ving - 1 cre | (10) 
TR 1 





Exaninin; th^ ۶ ^^ ir 


: TAS. aud 9, let 
Cie ctitical damping in the x direction 


LM m &, 


critical damping iu the y direction - 24m -e 
Substituting, Eq. 8 becomes; 


Cae = 


S he inf =e Oo 


(11) 
and since: 
Coc =2NMA = 24m pe 
Ri 7 
mk, = Sc, Le 
2 Y ka 
Eq. 9 becomes: 
7 e C7 i yo چاو دی‎ .. 
Re i t 7 | = (o) Au, - 6 n 8 (12) 
i ۱ 


Taking the damping ratios in the x and y directions to be equal 
Ge., ei/e, = ۰/6 it is apparent from an examination of Eqs. 10, 
11 and 12 that four independent dimensionless parameters define the 
problem, namely: 





_ me” 
6۱ I 

23 ے‎ Mm. 
FO 7 CIA 
un. be 


K 


( جح 2) < ( 2) = (4)م 


where k, is equal to or less than k,. 


Performing the above substitutions, 
Eq. 10 becomes: 


CG E E ec) ک‎ pra O -Y wo e| (13) 





.. e a Be , 

fongxA) f ee (6) ee 6 46 سم‎ (14) 
and Eq. 12 becomes: 
y + P(4) 4] PG) 7 n PG) y = (8) ow e - ۵ و همم‎ (15) 


Equations 13, 14 and 15 are the dimensionless equations of motion 
which will now be arranged in suitable form for solution by numerical 


methods. 





3. Numerical Solution. 

The Runge-Kutta fourth order method of numerical integration was 
selected to solve the equations of motion. This choice was made on the 
basis of several considerations, namely: it is applicable to non-linear 
differential equations; it is "self-staring", i.e., only the functional 
values at a single previous point are required to obtain the functional 
values ahead; it is relatively simple to program; finally, since it is 
essentially the Taylor series solution through terms of order Ga Lt 
offers an acceptable degree of accuracy for the solution of the problem. 
Before proceding, however, two additional parameters will be defined in 
order to provide the equations of motion with greater flexibility and to 
reduce the computation time, | 

The torque parameter P(2) is a constant for most of the presenta- 
tion. However, to provide for the case where it is a function of speed, 
let P(2) be replaced by P(2)P(5), where P(5) is either unity for the 
constant torque cases, or a function of speed for the variable torque 


case. 


Also, let: 
re JS 
si Az 
P(t) = PCS? = 1 5 
ZN, 


In order to employ the Runge-Kutta method of numerical integration, 
it is necessary to reduce the three equations of motion to a system of 


six first order differential equations of the form: 


dui ١ 
Oy 


$ 
IX 


Let; 
O = (16) 


T (18) 





Substituting these variables and the parameters P(2)P(5) and P(6) 


in the equations of motion, Eq. 13 becomes: 


x= vr) FU) و‎ Pu) J! سيم‎ ð- 2 C^ 4. ol (19) 
Eql 14 becomes: | 

| 4 P4) - f - ut08 Lei سم‎ © )20( 
and Eq. 15 becomes: 

È + PDP B PDY =x pua = Oho 6 (21) 


Equations 16 through 21 constitute the system of first order differ- 
ential equations which must be reduced to the form: 
d L ۱ PE e d L 
qe =$( 4%, O, fal); ; Hu (x) ) 


Appendix I contains this rearrangement. The resulting equations are: 


)16( کے 
T SE (17)‏ 
(18) اجوہ 
M | PG) Pls) + PO) PC) PL) Pcs OLAN P(3) y un O z (225‏ 


PU) P2 pue - PU) FT am e ۳ an 


A =| - PU)a? ou. B cu^ 6 + PUD PRAY pin © سد ©ممع‎ (23) 
PU) J a O 28 +h* om 8 ~ P(A) P(6) B - O = 


PO Pls) te © - Pl dt om? © + PL) F(4) P(e) 6 pu @ + 


P) F3) y % o |/ | PO] 





=|? foe O — F(4) y = Kf — POyYa? av^ 0 gea + (24) 


PU) P(4) P(e) 6 sin O loo O + RAN) PS? sm O uO- 
PUJA cea O T RNP NO O Lt PaO + 
P(2) Plz) Au 1/1] | — P() | 


Having selected the integration method and reduced the equations 


of motion to the proper form for solution by this method, the program 


itself will now be considered. 








4, Computer Progran Design and Use, 

Several factors governed the form of the computer program, In the 
material that follows, the more important of these are discussed in de- 
tail, while those of lesser importance are mentioned briefly. Details 
can be found in the block diagram or in the program itself, Appendices 
III and IV, respectively. 

Failure to accelerate through the critical speed region can re- 
sult from either of two closely related conditions, insufficient driving 
torque orexcessive unbalance of the rotor. An insufficient driving torque 
results in a condition wherein the angular velocity of the rotor remains 
in the vicinity of that corresponding to its natural frequency of vibration 
for a long period of time. In this case the inertia forces, though initially 
small, eventually produce large lateral deflections which in turn develop 
counter torques opposing the driving torque. An excessive unbalance pro- 
duces large inertia forces as a result of the greater eccentricity of the 
mass center. The resulting lateral defiections also produce counter 
torques opposing the applied torque. It follows, then, that the success 
or failure of an attempted acceleration through the critical speed region 
is dependent upon the applied torque and the eccentricity of the mass 
center, other physical quantities remaining constant. Hence a plot of the 
eccentricity of the mass center versus the applied torque for various accelera- 
tion attempts will indicate points of either successful or unsuccessful 
accelerations through the critical speed regions. A boundary line may 
then be drawn to separate those points representing successful accelerations 
from those points representing unsuccessful accelerations. 

For a given system, the parameter P(1) is directly proportional to the 
square of the eccentricity of the mass center, while parameter P(2) is 
directly proportional to the applied torque. These parameters, the eccentri- 


city parameter P(1) and the torque parameter P(2), are the ordinate and 
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abscissa, respectively, of the curve which is the boundary between the 
areas of successful and unsuccessful accelerations. The boundary it- 
self is defined by pairs of points. Each pair of points represents two 
acceleration attempts, one successful, the other unsuccessful, differ- 
ing by five per cent of the torque parameter P(2), all other parameters 
being constant. It is now appropriate to define a successful and an un- 
successful acceleration. 

Referring to Eqs. 20 and 21, one of the natural circular frequencies 


of vibration of the dimensionless system is unity and the other is the 
Ra 
Ri 


latter is unity. On this basis, a successful acceleration was defined 


square root of the stiffness ratio, . The maximum value of the 
as an acceleration in which a (dimensionless) angular velocity of two 
(twice that corresponding to the higher natural circular frequency of 
vibration) was reached. That is, if a dimensionless angular velocity of 
two was reached, the acceleration through the critical speed region was 
considered successful and the program terminated. 

The possibility existed, however, that during the fluctuations in 
the angular velocity of an otherwise unsuccessful acceleration a speed of 
two might have been momentarily attained, upon which the program would 
have terminated and the acceleration would have been erroneously termed 
successful. To investigate this possibility, several random successful 
and unsuccessful accelerations were initiated, and graphs of dimensionless 
angular velocity versus dimensionless time were made, Two such curves can 
be seen in Fig. 3. For the successful acceleration the torque parameter is 
five per cent larger than for the unsuccessful acceleration, all other para- 
meters being unchanged. It is apparent that the maximum angular velocity 
for the unsuccessful acceleration is well below that value (2.0) which would 


terminate the program and falsely indicate a successful acceleration. 
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An unsuccessful acceleration was defined essentially in terms of a 
successful acceleration. Starting from a known successful acceleration 
the torque parameter P(2) was decreased in five per cent increments, while 
the time to reach a dimensionless angular velocity of two was recorded. 
In regions far removed from the boundary between the successful and un- 
successful acceleration regions, the time required for succeeding success- 
ful accelerations increased by five per cent. Im the vicinity of the 
boundary, however, the additional time required for a successful accelera- 
tion exceeded the expected five per cent, but in no case did it exceed 
thirty per cent. On this basis, the time allowed for each successive 
acceleration was twice that required by the previous successful acceleration. 
If a dimensionless angular velocity of two was not achieved within this 
period of time, the attempt was considered unsuccessful and the program 
terminated. Figure 3 also illustrates this point. The upper curve re- 
presents a successful acceleration in which a dimensionless angular velocity 
of two was reached in time t. The lower curve represents an unsuccessful 
acceleration ín which a dimensionless angular velocity of two was not 
achieved in time 2t. The torque parameters for these attempts differ by 
five per cent, all other parameters being constant. Hence these two 
attempts correspond to one of the several pairs of aquta sun nl, define 
the boundary between regions of successful and unsuccessful accelerations. 

For the undamped case, an energy balance affords a check on the valid- 
ity of the results. At any given time the total energy supplied to the 
system must equal the sum of the kinetic energies of rotation and trans- 


lation, and the potential energy stored in the springs, or: 


dO. d Yon? A 4 - 
= Au —— d4 — E e 25 
| 4 = | - 2 
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An equivalent statement in dimensionless form is: 


| <)۱( | (o) 


eo) e= 7 3 | He 
SS [Si ی‎ 


(26) 


For the undamped runs, this check was performed at the end of each 
acceleration; i.e., employing the terminal values of the quantities 
in Eq. 26. 

The choice of a suitable time increment was made on the basis of 
a comparison of the results from several solutions of the same problem 
employing various time increments. Details are given in Appendix II. 

A (dimensionless) time increment of 0.10 was selected. 

Since the objectives of this investigation included a determination 
of the maximum amplitudes of vibration upon acceleration through the 
critical speed region, this quantity, the resultant of the deflections in 
the x and y directions, was computed each time increment and compared with 
a stored maximum which had been determined earlier, the larger being re- 
tained. 

Provision was also made to terminate the solution at any point and to 
resume the solution at this same point at a later date. This innovation 
provided a considerable savings in computer time since any time period, 
however small, could be completely utilized, and no solutions were lost 
because of unexpected requests to release the computer. The termination 
and resumption was executed by writing the values of the (dimensionless) 
deflections and velocities in the x and y directions, the angular velocity 
and displacement and elapsed time on magnetic tape in binary form, and 
reading these same values back into the program by the setting of certain 
selective jump switches. 

By way of a check of the computer, the output tape unit and the tape 
itself, the same short run was performed at the beginning of each period 
of computer use. 
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The program output included the parameters employed, the results of 
the energy balance, the maximum angular velocity and the time of its 
occurrence, ind the final values of the deflections and velocities in the 
x and y directions, angular velocity, revolutions, and elapsed time. Also 
recorded were the maximum amplitudes of vibration in the x direction, y 
direction and the maximum resultant amplitude. For each of the three 
maximum amplitudes of vibration, the time of its occurrence, angular velo- 
city and the x direction, y direction and resultant amplitudes of vibra- 
tion were also recorded. 

In general, the investigation proceeded as follows. Selecting an 
appropriate set of values for parameters P(1), P(3) and P(4), the torque 
parameter P(2) was varied to define the boundary between the regions of 
successful and unsuccessful accelerations. Having defined the boundary, P(2) 
was then increased through an appropriate range of values to obtain the 
corresponding maximum amplitudes of vibration. This procedure was then 
repeated for different values of P(1), P(3) and P(4). 


The following section presents the results of the investigation. 
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5, Discussion of Results. 

The results of the investigation are presented as follows. 

First, the boundary between the regions of successful and unsuccessful 
acceleration is discussed for stiffness ratios of unity, 0.50, 0.75 and 
0.25. Then the curves of the maximum dimensionless amplitudes of vibra- 
tion are presented. Finally, the effects of added damping are considered, 

A rotor having equal spring constants in the x and y directions (P(3) 
= 1.0) and no damping (P(4) = 0) was investigated first. Figure4* indicates 
the location of the boundary between the regions of successful and un- 
successful accelerations as a function of the eccentricity parameter P(l) 
and the torque parameter P(2). This curve is essentially a reference for 
the boundary curves that follow since unequal stiffness and added damping 
are more conveniently discussed in terms of their effects on this particu- 
lar result. 

Figure 4 shows that at successively smaller eccentricities, smaller 
applied torques are required for successful accelerations, since the 
smaller resultant inertia forces produce smaller counter torques. 

Decreasing the stiffness ratio to 0.50, the undamped boundary curve 
was found to be as shown in Fig. 5. For any non-unity stiffness ratio, 
the system will have two critical speed regions; in this particular case 
one is in the vicinity of a (dimensionless) angular velocity of 0.50 
and the other is in the vicinity of unity. It follows, then, that an 
unsuccessful acceleration can result from an inability to accelerate 
through either the lower critical speed region or the higher critical 
speed region. Figure 5 illustrates both of these conditions. The bound- 
ary between the regions of successful and unsuccessful accelerations is 
seen to vary between two parallel limits. In the vicinity of the higher 
limit, the unsuccessful acceleration points defining the boundary represent 
failure to accelerate through the higher critical speed region. Similarly, 


*Figure 4 can be found on page 26. 





those unsuccessful accelerations defining the boundary in the vicinity of 
the lower limit represent failure to accelerate through the lower critical 
speed region. 

The limits mentioned above have a definite relationship to each other, 
For a given eccentricity, the ratio of the lower limit to the upper limit 


is very nearly equal to the stiffness ratio, kj/k,. This can be explained 


M m 
TRI 


It is apparent from this dimensionless grouping that a system which has a 





by examining the dimensionless torque parameter P(2) = 


spring constant of, say, Ak, will require a torque of AM for a given 


1 
eccentricity. In the vicinity of the lower critical speed region of the 
system under consideration, the system behaves essentially as a single de- 
gree of freedom system having a stiffness k, = 0.50k,. Hence the torque 
required for acceleration is 0.50M. Since the boundary curves are plotted 
versus P(2), which is a function of Ki: the substitution of 0,50M for M 
results in a decrease of the torque parameter from P(2) to 0.50P(2). 

The fact that a portion of the boundary curve is in the vicinity of 
the lower limit shows that, in this region, a successful transit of the 
lower critical speed region will result in a successful transit of the 
higher critical speed region also. This condition may be due to a 
favorable phase relationship between those vibrations induced at the lower 
and higher critical speeds. That is, those vibrations remaining in the 
system as a result of the transit of the lower critical speed region exert 
forward torque opposing that developed by the vibrations induced in the 
higher critical speed region to the extent that a smaller torque parameter 
suffices for a successful acceleration. 

For a stiffness ratio of 0.75, the undamped boundary curve was found 
to be as shown in Fig. 6. The general shape of the curve is comparable to 


that of the stiffness ratio previously discussed, 0.50; however, some dif- 


ferences do exist. 
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The boundary is seen to vary between two limits much the same as 
the boundary of the 0.50 stiffness ratio case; however, for a given eccentri- 
city, the ratio of the torque at the lower limit to that at the upper limit 
is influenced considerably by the close proximity of the two natural fre- 
quencies of vibration. In the higher torque کے‎ this influence is 
quite apparent, and as a result, the limit ratio differs somewhat from the 
stiffness ratio. In the lower torque regions, agreement between the stiff- 
ness ratio and the limit ratio is quite good since there are many revolu- 
tions of the rotor separating the critical speed regions and their inter- 
action is reduced considerably. 

For a stiffness ratio of 0.25, the critical speeds of 0.50 and unity 
are sufficiently displaced from each other that little or no interaction 
of vibrations is experienced. The counter torques developed in the higher 
critical speed region are larger than those induced in the lower critical 
speed region due to the larger angular velocity. It follows that an un- 
successful acceleration represents a failure to accelerate through the 
higher critical speed region. The unsuccessful acceleration points defining 
the boundary curve shown in Fig. 7 represent just such a failure. 

By way of a brief summary, the undamped boundary curves described 
above, i.e., for stiffness ratios of 1.00, 0.50, 0.75 and 0.25, are 
plotted on a single graph, Fig. 8. 

Except as noted below, the starting phase 9 = 0 was used. When the 
lateral stiffnesses are equal (P(3) = 1), the resulting circular symmetry 
assures that neither the location of the boundary nor the maximum vibra- 
tion amplitude will be affected by starting phase. Even with unequal 
stiffness ratios, it is reasonable to expect that the effect of starting 
phase is negligible if the rotor completes a large number of revolutions 
before reaching the lower critical speed. This condition is met for 


small values of the torque parameter. 
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An anomalous behavior in the higher torque region did result ina 
limited study of the effect of starting phase. With a stiffness ratio 
of 0.50, it was found at P(1) = رت‎ that the minimum value of )2( 
for successful acceleration was 46.0 per cent less than at immediately 
adjacent higher and lower values of P(l). Further runs were made at 


P(1) = 2.3010 > 


with starting values of 0 - 45, 90 and 135 degrees and 

it was found that the minimum P(2) for successful acceleration approached 
the results obtained for neighboring values of P(1). Because the principal 
engineering interest in these results is in selecting a combination of 
parameters that will assure a successful acceleration, regardless of start- 
ing phase, the "stray" point at P(l) = 230x100 obtained with O - O was 
ignored in plotting the boundary curve of Fig. 5. 

Figures 9 through 12 are the undamped maximum amplitude curves for 
the stiffness ratios investigated. These curves indicate that, for a 
given eccentricity, successively larger applied torques result in de- 
creasing amplitudes of vibration. These smaller amplitudes are to be ex- 
pected since the larger torques result in a shorter period of time in the 
critical speed regions. 

These curves further show that in regions of successful accelerations 
removed from the boundary between the regions of successful and unsuccessful 
accelerations, the maximum dimensionless amplitudes of vibration are in- 
dependent of the eccentricity parameter P(l); or, for a given torque para- 
meter P(2): 

R/e - Constant 
from which: 

R = (e) (Constant) 
Therefore, the maximum amplitude of vibration is directly proportional to 
the eccentricity e of the mass center. This is to be expected since these 
vibrations are a consequence of the inertia forces which, in turn are direct- 


ly proportional to the eccentricity of the mass center. 
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The unusual variations in the amplitude curves of Figs 10 and 11 
are a result of the transition of the boundary curve from one of its 
limits to the other limit. 

Figures 13 through 20 indicate the effects of the addition of one 
per cent of critical damping on the boundary and amplitude curves, With 
one exception, shown in Fig. 15, these curves reveal a decrease in the 
maximum amplitudes and a decrease in the torque required for a successful 
acceleration. The decreased amplitudes result in smaller counter torques, 
hence acceleration through the crítical speed region is achieved with a 
smaller applied torque. The larger shift of the boundary in the lower 
torque regions is a result of the larger number of cycles during which 
the damping has acted. 

Figure 15, however, shows that the effects of added damping are un- 
favorable for a considerable range of values of eccentricity and torque 
parameters. It is possible that the added damping, though reducing the 
individual amplitudes of vibration induced in the lower and higher critical 
speed regions, could have also resulted in a particularly unfavorable phase 
relationship between these vibrations, As a consequence the resulting 
amplitudes and counter torques could have been particularly large and a 
larger applied torque would be required for successful acceleration. 

As was the case with Figs. 10 and 11, the unusual variations of the 
amplitude curves of Figs, 18 and 19 are a result of the transition of the 


boundary curve from one of its limits to the other. 
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6. Sample Prcblem, 

The computer program is not limited to accelerations with constant 
applied torques. The following is a sample problem illustrating the ex- 
tension of the solution to a system having a variable, rather than a 
constant, applied torque. 

An unbalanced rotor, its supporting shaft and drive sheave have a 
combined mass m = 12.0 lbm, and a moment of inertia about the mass center 
I, = 48.0 Ian The mass center of the system is displaced a distance 
e = 0.02 in.from the center of rotation. The stiffnesses of the shaft and 
its supports are the same in the x and y directions and are ki = k, = 100 
lb£f/in. The drive motor for this system has a starting torque M = 40.0 lbf- 
in- which decreases linearly with angular velocity to a value of zero at 
55.2 rpm, or twice the critical speed of the system. If damping is negli- 
gible, determine the maximum amplitude of vibration. 

The axis of the moment of inertia must be translated from the mass 
center to the center of rotation. Since: 


t= le Mer = 


T= (480 Ibm - in?) + (12.0 Ibm)( 0.02 in) = 49,۵ ٣ 


Defining the parameters: 
2 


Dreamer 


1 
AO کر‎ = 1.00 x one 
(48.0 Ibm in?) 
Similarly: 
Mm 
FE) m em 
(2) IR, uu 
2) = (40.0 lbí-w.X 12.0 Ibm 0.100 


(48,0 lom- i4?Y 100 lbf /in) 


This value of torque decreases linearly to zero at a (dimensionless) angular 

velocity of 2.00 (twice er speed), therefore the modifying parameter 

P(5), which has been unity for the constant torque case, now becomes: 
Pe)=(1- —) 


and the dimensionless torque terms in the equations of motion become: 
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Also: 

Pi mee ee. وت‎ 

(3) E OQ 
and, since damping is negligible: 

0 ٤ 
Inserting these parameters into the computer program, the resulting solution 
indicates: 

R 5 1.49 

= 
from which: 


pec YO im. 


where R is the maximum resultant amplitude of vibration. 
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case). 
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of the boundary between regíons of successful and unsuccess- 
ful accelerations, stiffness ratio P(3) = 1.00, P(4) = 0.02. 
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APPENDIX 7 
CONVERSION OF EQUATIONS FOR NUMERICAL INTEGRATION 
The system of first order equations to be rearranged is restated 


for convenience. 


(16) 


1 an 


fen (18) 
A= POJP(5) + PU) | itas 6 fee e | (19) 
a fo DITS uuu EE ee (20) 
6 + P(4) P(6) A + Hamm: ےت‎ Low 9 (21) 


Equations 16, 17 and 18 have the required form: 


UE ( 
HE = fl, O) O) ره‎ Am) 
m f 2 ^ ) 02 ) ) 7 
It remains, then,to reduce Eqs. 19, 20 and 21 to this same form. 


Rearranging these equations: 
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Employing Cramer's Rule, 
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APPENDIM TI 
WE Paar ION OF TIME INCREMENI 
The selection of thé tine increment h was made by selecting several 

representative sets of parameters and obtaining the corresponding solu- 

tions for various h. This afforded two comparisons, one of the energy 
balances and another of the resultant deflections, Then the torque para- 
meter was varied until the boundary between the regions of successful and 
unsuccessful accelerations was defined. This afforded à third check,  Re- 
presentative results for one set of parameters are presented in the accompany- 


ing table, The time increment of 0,1 was selected as the reference since 


| 
smaller time increments produced no chänge in any of the three comparisons 
made for each set of parameters. 
Per cent difference as compared 
with results for h = 0,1 
EE II ` ات‎ 7 


۱ 
. 00 0.00 5.00 


O 
N 
O 


po 0.61 0.00 0.00 
0.4 223 2.0 0,00 
no 1,82 7.20 0,00 
Effect of time increment h on resultant 
deflection, energy balance and location 


of the boundary between regions of 
successful and unsuccessful accelerations. 


47 





Start| 


[sas ےت‎ 


E 1 
ior test rw | 





wu s] 
= Set, 261 lal va ues - 


n Ñ. — Ñ - nm I‏ — ' ا — امم یلا 





- فاگ سس‎ 
| ;valuate integrals 


Compute resultant 
deflection R 













Replace with 


272۳7 و‎ 
store R. ٦ 














| Replace with 


۱ و larger Y‏ 
۵ و 5 و و store‏ 





15  ) ٦ 







Replace with 
larger R, . 
store n,1,5,9 








(p/P, snared 


PINE SER 
outputs 





Set new 


parameters 


© 








AD DES 


COMPUTER PRU. RAM 


7 
مہ‎ 
> 
~ 
کا‎ 
M 
e 
~ 
N 
— ren 
کت‎ a ےی‎ 
> pe] 
— m ol) 
o سم سے‎ XX om 
TE — een ہو سے‎ 
— >= — — N) 
ھ8‎ ^t ليك ہے‎ = 
ml) سے‎ en 
^ NA Hn مس‎ 
— — سه‎ ta > N) 
© ند‎ ED — — 
BS a de CN Bullet 
اھ‎ > ~ =< (nun 
eg € mu nn Zu 
© — © ka Came 
> — < —— 8ھ‎ 71 
E سم‎ e سم‎ PO NSN O ہم‎ 
COSE — سح‎ vie CC e o ° ée >> مه سمہ سے مہ ےم — — — سپ سد‎ 
= — DOCH پچ وا‎ ZF 0 $14. 1.1400 EE L et eg EE >> >< >= 
LANZ » A SS RAN O O IA A ۱! ۱۱ H mH NK NAN HE H ۱۱ 
کم‎ ig TD دی‎ HW o H J! LIKIDE jH H M ill ——— “cC OQ ><: > ED CO A e رس بح جح‎ mn O es 
a La) som ne TE جم‎ o A BR LE AA AA AO qup 
ے ہے ے بب‎ ب٣‎ ہ٣‎ E wl WS LI wl E Léi Pr CN) TLIO O ت‎ iwi IXN ON Dun uuuuuu =— 
G Kd CD a || سب سب ست سي‎ IU LU IL De. Ze un || LUN I HLA لا لا‎ „IDOL UON 


سج U XI NL mm > > > > > LD‏ لچ سر سح کچ > > > > > > <> ہی LEE G D‏ کت کہ ات نم ام بے لد کے هل نا D.‏ 


5 © Ec - 
Le e سم‎ KN == ریہ‎ ze Sc? 








91,7? 


2.) 


3 
pon 
پر‎ 
VY) 
——  — oO 
N SO =N D للا‎ 
— —— صد‎ (€) 


u 
> C > >< <N > > x > > Z کیہ‎ 
! || fe || زلما‎ kee Zei e 
= LI — مط کہ > سا لا کے‎ o لاا کے ا‎ 
رح‎ AB UO Re mM OUO zum sr Y IOS mR e LUS wi EC E 
DZ Zum Zee UNI ےی‎ — O — 
سی مھ سے ( ےہ اکر امت لع ذم‎ SL تم تلا مه بالط‎ es E 
>>> > >>> > > — Q. t ZL T. L Zt. u 2 سنا © سلا‎ 2 


H H |( tl‏ !ا 


w w a ` ےہ ہد صےد‎ w تسس‎ 


NË 


M 


Mm 


>< 

ra 
نے‎ 
O دم‎ 

e e. 
ںہ سم‎ 

H N 
Loan] 

Om سے‎ 
na — 
ہم‎ 
E 
w en 

o — 
¡a ۰ ید‎ 
SI 


NUN وق و قل" ور ا‎ I 
Le ON reg bes K سس‎ (Doc 


N 


),8,8 


">" > 
N و‎ 


2 


ANS 
INN 


O re 
EU puc EE 
سے سم ےم‎ Ni; 


C Sn 


IN 


=> 
-7 
ee 
SS Saj 

~ 
N سب‎ 
9 a 
= N لن‎ 
— S CH 
د‎ x جج‎ 
بر(‎ *— — 153) 
مت سس‎ N Se ex 
كد‎ > ER 
On 6ت‎ whe سے‎ 
> سوب‎ x ےج‎ 
SEEN, xt N 
CN — > CN: N en: 9 
< aeu Se 
سج ت بے‎ ~ E 
r rm Dën LI e >e LL. 
Li +— a 
Dm ALIN PN) T سے‎ 
ہے ہے‎ S S || x nO t 
“ e eh سد ي‎ > 
<= StL N 
لہ‎ UND oN ١ 
۱ CG ai ~u = 
Lee ات‎ > U, = 
D کک سی و‎ EE Le 
A ISR > 
ای سك‎ 60 ۱۷۷۷۸ GG 
XKX I سح‎ o aem للا‎ 
اہ کے‎ m کی سو چا تہج‎ 
Don SK mm re اس مد‎ Nu 
xD ud. f C. x | “NCNM 
EN muss IR ER AX e YO 
ODT mw mm DIL ® 


`C —! =‏ پچ AH‏ < صہ لے 
SU ü CN mm — — — an‏ عدم 


jet» 


LL — — ید‎ — — > `” Ə — >x G. + —— — =. 


O = ذه‎ 


SO Oe ee 


em OD) m NT ZD> سم‎ 


> O. 
مس لیا‎ x ل‎ 


Ge EE Ee 


— 


Zu T CN OS I] >= a 
ec oeil 
vd CN mme 1ا‎ 


Ni Hat UN C >< 


verad 
یں‎ AN f 


eS wp + 1‏ ہس 


1l, 85 


1) 


n س‎ lt ies N 
ev lf ` An £ C = (| TI حل‎ 


Erna rt WV b 
—4 — کک“‎ — — 


N TE ور و‎ i مہہ ود ہہ‎ Z > > AAA 
F —— ملاح مد رد‎ eet A 
e LN ES, D وي‎ X Zë ee LL Ye 
<A „In ےھ‎ SNN Meow ht «mec mn 


= 


— 


= 


حم 
أيهم 


U^ 
Pa 


bes kees IL zm H — Q.‏ ہم 
سای سر ےم سج نے جج 
ÑZ e‏ ری سس 


can P 


cu 
TINTE 
DEN سس‎ Gär 
Se — 
= m 


نے < 
تلا نلا با 


P الى را ا‎ y; pl) 


NO II N OO n UJ. 


ب 
on‏ 
Ta‏ 
emm‏ 
a‏ — 
s N‏ — 
سے سم “Oy‏ 
زج مس و — 
el E‏ 
نل لان ںیم 
Ze 20‏ 
ZO} wo‏ 


OA OAC — y A — 
>=> AA a AED SAA 


— N] 


ہی LL St LUC Co ee‏ لالہ © سم كه 


rat Kt D DELL ا ساي‎ 


KS 
DICH em 


مم |" 
~ 


TS 
a 





D 
ES = Së 
E om — 
—— — سے‎ vu Lë — GE 
DO — ~a” ~ پر ہہ‎ ۰ 
— = >Ə NO coe Oy —— < 
> CO TT ~ ۰ 
x ہم‎ x O. x جح‎ ES 
— OM i = ~ 
T E m nO po E 
~ — N ہم‎ m 
CL ی‎ x بل‎ e + 
E A یلم عیب‎ zt RA سے‎ 
oi Qe — N^, SZ 
سے‎ ce y — (2) iS ch 
= + NO ET > e 
E (Ox سے سه‎ e => 
ws لن‎ — ~ 
جم‎ > N یه‎ 
a OO س‎ > = 
NN | = x سے کت‎ a 
E ~ Na — e e 
xoc DUE COO e LJ N 
us Mo دی کے حا ری‎ Ws — 
کے تا‎ — — — rns — (257 ~ 
zT Se E E aS p T CN 
SEH — a. — مس‎ j P. — T > = 
مد‎ N~ >~ Q(Q > > m x 
zf kk + X >= Ee < 
— ^o — — — ar | x 
بس — = — => سح و نب‎ ° 
yx — ine ہا مہ سی‎ (7T e KE س‎ GN: 
مه‎ A سے سے‎ — a. اک‎ ^ m >< xz SR 
en pd Mi Ne SJ | ب‎ Z ^ سے‎ - 
ووس يي سم ماس‎ | SO eo u= م‎ — 
>x — vn as mas ہے‎ cc | _ 
aX x سم سے سس‎ x — o ہم < ل اكد جح‎ Y 
> =< we oe ee No ہے‎ Lo => A ہے جح وس‎ 
سح یج | ره ہے‎ emo Z= > üJ (Cy c سب‎ 
ao ہم = مس — = اه سح‎ Xx 7 + 
ches pou — تب دج ها و‎ 
> >- سح سب‎ X + * — ck e eu were UO >» = 
ح ات‎ > —=— ez ep مس‎ Léi NA eee See — 
ES om M IA LH ui 2 “SO OOD et UY 
۱۳ 0 — ناسح — ہج ات‎ NC S 2^2 x ۱۵۵ ee ہف‎ 
سک‎ — —X aa ut Pae Ç > (e i Y > A a rd UN ¡e A 
a mmc و سس‎ s = cu gcc UM 
u تن ہے یہ ہت‎ Sa ice لب‎ > Z = =y 2 Z <( >— Z — Z > 
= ER >= >= >> — — — CN! IE ہے‎ N عت‎ æ ' mat K kaa al) e! 
اانه "اسمن عن ص ره اس رح 2 بت‎ u u and N nein 


LL LL en HD O AmQ‏ سر مس 
ب > > -> DNT ZNT z‏ 


UN uu ۱ a>‏ ا بن 
I sn‏ ہے س وہ وہ تّ DIS‏ 

ص > ^ HY Ao Im H QO Z= H H iH‏ زا اا لاس بے ری کے كك 

E ہم حم وى سر زی تم‎ dE ہم‎ | Een Y DAA ern یس تی ید ہچ جح )ہم کو یں‎ SD 
pr = ft PENA ANOS سم‎ SS LN) س ص چ‎ n AO CN < OC Z NAY LI مب‎ j Ze Da 
= مہم رے‎ NM Zum r Vu A Deut) ١ے”‎ ux = — — 0 >x — C) C U < C; >x G > لا‎ 
ات مت ف ف رك رك رك هم حي انم‎ OCO + حت ع لب‎ SIO OO Sa >< l S D p a pusi 
' Deren rend 


سسا ص 


مو — 


)جح 


S N مک الچ‎ 


un 
شا‎ 








L 


POINTS 


b 
| ?2z SSI09RS | ?27ss322n$Nn|  c04ss220nsS | PI SIREN) (0d Í 


m e 1 
7 






(g2¢Wvd ) Z0'0 =(t)d | (eerwyenn) 0900 =(v)d 
ST Q = (€) POF (Id FO SIENTEN 





E ہے‎ 








1-9 * eb'b | 4-0 x Tbb Q£ 0۷7 | ۲0 7 1 
AO. A A NE vo | 
_ و‎ 2٤۶۶ he$| g-0 RS | 0m 7 | saru 00% | 

1 








G 00 ۶ (۵ 
001 = (Ee)d 


) 020۳7۷۷7۸۷ ( 2000 2 ۵ 
WOT (Gd 30 S o m 











Points Defining the Boundai 


_ و‎ mel ar جع‎ 
| z- 9 * 7 ٤۷۷ ) Let gel ۶'٥٦٣ 
1 77755922757 ١/7505 ١ 7 


(Bumva)rToo= (WI | (azdwvenn) 00'0 = (P)d 


FLO =2(E)F 804 (¿a (۸ 





| QO F | 
OC 
200. 

_O/ ۸۲ ںی‎ ! 
2 

p-r 00'E t 
00 | 
40,۱۱ 
5 N 
Op eT 
272۶ 
OS" C 
OB b | 
C 


> 


— 


OO: 
(712 


- m 
Q 
x 
9 
O 
Ga 





Cont, ) 


undary 


ths Y 
cn D 


ing 


Points Defir 










eee cz on SCHT ven. 
m 
















APPENDIX V 


MAXIMUM 7ھ‎ (PLITUDES OF VIBRATION FOR SUCCESSFUL 
ACCELERATIONS THROUGH THE CRI "TICAL SPEED REGIONS 
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